Level-crossing spectroscopy of the 7, 9, and IOD5/2 states of ^^^Cs and validation of 
relativistic many-body calculations of the polarizabilities and hyperfine constants 
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We present an experimental and theoretical investigation of the polarizabilities and hyperfine 
constants of Dj states in ^^^Cs for J = 3/2 and J — 5/2. New experimental values for the hyperfine 
constant A are obtained from level-crossing signals of the (7,9,10)D5/2 states of ^'^''Cs and precise 
calculations of the tensor polarizabilities a2- The results of relativistic many-body calculations for 
scalar and tensor polarizabilities of the (5-10)D3/2 and (5-10)-D5/2 states are presented and compared 
with measured values from the literature. Calculated values of the hyperfine constants A for these 
states are also presented and checked for consistency with experimental values. 

PACS numbers: 32.10.Dk,32.10.Fn,31.15.Ar,31.25.Jf 
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I. INTRODUCTION 



Level-crossing spectroscopy in an electric field has been 
shown to be a useful technique to determine atomic prop- 
erties. Already the first experimental studies of reso- 
nant signals at pure electric field crossings of magnetic 
imp components of certain hyperfine (hfs) atomic lev- 
els F at nonzero electric field [ll, [3, 0| and their further 
development by applying two-step laser excitation Q 
demonstrated how this technique could be used to ob- 
tain atomic properties. The method makes use of the 
fact that the electric field values at which magnetic sub- 
levels niF cross in an electric field depend on the tensor 
polarizability a2 and on the hfs constants. When the 
electric field is scanned and laser induced fluorescence 
(LIF) of definite polarization is observed, these cross- 
ings are associated with resonance behavior in the LIF 
signals. When the separation between crossings is large 
compared to the widths of the resonance signals, as in the 
^L)3/2 states of the -'^^'^Cs atom (see Fig.[T^), they lead to 
rather well-pronounced resonances in the observed fluo- 
rescence. Moreover, these resonances correspond exactly 
to the level-crossing points under appropriate experimen- 
tal conditions. Such resonances were used to measure the 
tensor polarizabilities a2 in the 7,9D3/2 states of ^^^Cs 
atoms [J] , in which the magnetic dipole coupling hfs con- 
stant A had been previously measured with good pre- 
cision, and the electric quadrupole hfs constant B was 
assumed to be negligibly small [5] . 

Such measurements become more challenging, how- 
ever, in the case of the states of cesium, since 
there are many closely spaced crossing points of mag- 
netic ±mp components (see Fig. [TJb-d). As a result. 



the level-crossing signals overlap and no longer contain 
discernable resonances. In this case, reliable values for 
atomic properties can be extracted only by means of a 
very detailed and accurate theoretical description of the 
observed electric field dependence of the signals as a func- 
tion of atomic properties and experimental conditions. 
Such theoretical descriptions have been developed and 
tested in connection with the states [J. Neverthe- 

less, the level-crossing technique cannot be used at this 
time to improve the knowledge of the tensor polarizabil- 
ities a2 of the nD^ j2 states because the extant measure- 
ments of the hfs constant A contain uncertainties on the 
order of 30%. The small hyperfine interaction, especially 
for n> 7, makes them difficult to measure @, 0] ■ 

The first value of the hfs constant A of the nDc^/2 
states of ^'^^Cs was obtained with measurements of the 
widths of optical double resonance (ODR) signals in 
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the Paschen-Back region. The results for the 
and IOD5/2 states were -0.5(2) MHz and -0.4(2) MHz, 
respectively Q. These values were improved through 
level-crossing spectroscopy in magnetic fields, yielding 
-0.40(15) and -0.30(10) MHz %. The authors com- 
bined these data with previous ODR measurements Q 
and presented the weighted average as —0.45(10) and 
—0.35(10) MHz. They concluded that the quadrupole 
interaction can be completely ignored when fitting the 
experimental data. For the 7D5/2 state, Bulos et al. es- 
timated the A value to be -1.7(2) MHz from the ODR 
signal width Q . The drawback of the ODR experiments 
on the 7,9,10D5/2 signals is that it is necessary to use in- 
direct cascade transitions to observe the nDc^ /2 signals be- 
cause of the presence of scattered light at the nDs /2 — > 
6P3/2 fluorescence transition 0, Q 

The tensor polarizabilities a2 for the nD5/2 states in 
cesium under discussion are known with far greater preci- 
sion. For the IOD5/2 state, 02 has been measured by Xia 
and coworkers to a very high precision of about 0.3% at 
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6815(20) X 10^ a.u. [13]. For the 9D5/2 state the 02 value 
is measured with ca. 5% accuracy at 2650(140) x lO'^ 
a.u. by means of level crossing spectroscopy in combined 
electric and magnetic fields [a, For the 7D5/2 state, 
there exists a measured value of 129(4) x 10'^ a.u. pre- 
sented in Ref. [l3|, which, however, should be verified, 
because it differs significantly from the theoretical value 
of 140 X lO'^ a.u. [13]. Furthermore, a more recent mea- 
surement of oi2 for the 7D3/2 state [4] was closer to the 
theoretical estimate of 13] than the measurement of [T^. 

The situation with the electronic structure calculations 
is similar to the experimental situation. Rather good ac- 
curacy has been achieved for theoretical estimates of the 
tensor polarizability , as can be seen from the fact that 
the calculations of [ij] for ^'^^Cs agree with very accurate 
experimental data for the (10-13)D3/2.5/2 states [l^. De- 
spite this precision for the polarizability, the estimates of 
the hfs constants are poor and can hardly be evaluated 
reliably, for reasons that will be discussed below. There- 
fore, there is a need for more accurate values for the hfs 
constants of the states. 

In order to determine the hfs constants A from our 
measurements oi mp sublevel crossing signals in the 
7,9,10D5/2 states of cesium, we used the following ap- 
proach. We fit the measured signals with calculated 
curves derived from simulations, which had been devel- 
oped and tested in [3] • With the tensor polarizability a2 
fixed, these fits yielded the hfs constant A. To choose the 
proper value for a2 , we performed an all-order relativistic 
many-body calculation. 

Section II contains a description of the experiment, fol- 
lowed by a discussion of the simulations used to describe 
the measured signals. The all-order relativistic many- 
body calculations that provided the values for the tensor 
polarizabilities a2 are described in section III, and the 
values for a2 obtained from these calculations are com- 
pared with earlier experimentally measured values. In 
section IV we discuss the application of these calcula- 
tions to estimating the hfs constants and compare them 
to the results of previous experiments. In section V we 
show how to use our experimental results from section II 
and the calculated tensor polarizabilities from section III 
to estimate new values for the hfs constant A. 



II. EXPERIMENT AND DESCRIPTION OF 
SIGNALS 

A. Method 

The premise of level crossing spectroscopy is that the 
spatial intensity distribution and polarization of the laser 
induced fiuorescence produced when an atom is excited 
depends on the coherences between different magnetic 
sublevels mp of hyperfine levels F. Such coherences arc 
destroyed when the degeneracy between different sub- 
levels is broken in an electric field. However, in the case of 
linear polarization, they can be restored when sublevels 



with Amp = ±2 cross at certain electric field values. 
Figure [1] shows the hyperfine level-splitting diagram in 
an external electric field for the 7D3/2 and 7D5/2 states 
of cesium. This diagram is calculated by diagonalizing 
the Hamiltonian, which includes th e hy perfine and Stark 
interactions, in an uncoupled basis [14| . 

When applying the method of level crossing spec- 
troscopy to the study of the nD5/2 states of cesium, 
one encounters two difficulties not present in the case 
of the nD^/2 states. The first difficulty is that the ^05/2 
hyperfine manifold contains seven level crossings with 
Amp = ±2, whereas the nD^/2 manifold contains only 
two (see Fig. [T]). The large number of level crossings in 
the nDs /2 state wash out the sharp resonances that could 
be observed in the nD^/2 state. 

The second difficulty is that in the case of the ^05/2 
states, after the two-step excitation 6S1/2 — > 6P3/2 
— > f^D5/2 (see Fig. [2]), it is necessary to observe the 
fluorescence from the nDQ/2 — ^ 6P3/2 transition. Thus, 
scattered light from the exciting laser constitutes a high 
background that must be suppressed. Figure O shows the 
level excitation scheme. 



B. Experimental details 

We studied cesium vapor at room temperature in a 
glass cell. The experimental setup is essentially the same 
as in Ref. [!]• We could apply an electric field between 
two transparent electrodes inside the cell, which were sep- 
arated by a 2.5 mm gap. Figure [3] shows a schematic 
diagram of the experimental setup and geometry. The 
most crucial detail of the experiment is the relative ori- 
entation of the electric field and the polarization vectors 
of the linearly polarized laser radiation. The first laser, 
which excited the 6S1/2 — *■ 6P3/2 transition, was polar- 
ized with its polarization vector Ei parallel to the dc elec- 
tric field £, which was along the z-axis. The second laser, 
which excited the 6P3/2 — > 7,9,10D5/2 transition, was 
sent in a counter propagating direction and was polarized 
perpendicular to the first, with polarization vector E2 
parallel to the y-axis. We observed the laser induced flu- 
orescence (LIF) at the nD^/2 — > 6P3/2 transition along 
the z-axis through the transparent electrodes. A linear 
polarizer selected the intensities of the LIF polarization 
components along the x or y-axes Ix or ly . Since the LIF 
was observed at the same frequency at which the second 
laser was operating, it was necessary to suppress care- 
fully the scattered light by means of diaphragms. The 
scattered light accounted for between 30% and 50% of 
the measured signal. We checked that this background 
remained stable during the measurements and subtracted 
it from the signals. The LIF passed through an MDR-3 
monochromator with 2.6 nm/mm inverse dispersion and 
was recorded with a PMT in photon counting mode dur- 
ing one second time intervals. 

The first laser was always a diode laser (based on an 
LD-0850-lOOsm laser diode) and was tuned to excite the 
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FIG. 1: Hyperfine level-splitting diagram in an external elec- 
tric field for the (a) 7D3/2, (b) 7D5/2, (c) 9D5/2, and (d) 
IOD5/2 states of Cs. Circled points indicate level crossings 
with Amp ~ ±2. 




F=3 tilt 
FIG. 2: Level excitation scheme. 

6^P3/2 state from the F = 3 his component of the ground 
state. We chose to excite from the F = 3 level, because 
in this way we could avoid the F = 6 level of the nD^ /2 
final state. The F = 6 level contained no level crossings 
and thus would contribute only background. We took 
advantage of a sideband of the radiation of the first laser 
in order to achieve broadband excitation. 

For the second excitation step, we used a diode laser 
(based on a Hitachi HL6738MG laser diode) in the case of 
the 7D5/2 state and a CR699-21 ring dye laser with Rho- 
damine 6G dye in the case of the 9D5/2 and IOD5/2 states. 
The second laser was operated in single mode regime. 
We recorded data at different values of the detuning of 
the second laser in order to compare the results obtained 
at different detunings with simulations. A HighFinesse 
WS/6 wavemeter allowed us to measure changes in the 
lasers' detuning with a resolution of 30 MHz. However, 
in general we operated at the detuning that maximized 
the fluorescence signal. When the second laser was the 
diode laser, we jittered its output frequency over a range 
of approximately 1.2 GHz by applying a sawtooth wave 
with a frequency of tens of Hertz to a piezoelectric crys- 
tal mounted to its feedback grating. The laser power was 
of the order of a few mW, and the laser beam diameters 
were approximately 1 mm. 

The electric field produced in the cell was cali- 
brated with measurements of level-crossing signals for the 
IOD3/2 state of cesium as in The level-crossing res- 
onance positions obtained with our cell were compared 
with the crossing points calculated from the tensor polar- 
izability of [10] and the hfs constant A of [HI] . The overall 
uncertainty on the electric field magnitude was estimated 
to be about 1%. 



C. Experimental results 

We plot with markers the measured LIF intensity as 
a function of the electric field for the nDi^^2 states of ce- 
sium in Figures UHHl Signals for different experimental 
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FIG. 3: Experimental geometry. 



geometries are plotted. We label the experimental geom- 
etry as zyy or zyx^ where the first and the second letters, 
z and y, denote the orientation of the polarization of the 
first and second lasers, Ei||2: and E2II?/, and the third 
letter, x or y, denotes the polarization direction of the 
observed LIF. The solid line in the figure indicates the 
results of the simulations that are described below. As 
inputs to the theoretical model, we used the tensor po- 
larizabilities calculated with the relativistic many-body 
approach described in section III below. 
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FIG. 4: Experimental results for lT)^i2- 



As can be seen from Figures SI [SI and [5] there are no 
well-defined level crossing resonances. Nevertheless, a 
curve with multiple features is obtained, and these fea- 
tures can be fitted with the results of a simulation based 
on a theoretical model. This simulation is described in 
the following subsection. The fit involves adjusting the 
hyperfine constant A and those experimental parameters 
that we could not measure absolutely, such as the laser 
detuning. We fix the tensor polarizability at the val- 
ues that are obtained from the calculations described in 
section III. 
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FIG. 5: Experimental results for 9D5/2. 
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FIG. 6: Experimental results for lOD 
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D. Signal description 

Since well-defined resonances are no longer present in 
the signals of the D5/2 states, the data can be interpreted 
only by means of simulations based on a detailed model. 
Such a model was elaborated in detail and verified in 
a previous publication fJ|, so we will only outline the 
approach in what follows. 

The model describes atoms that interact simultane- 
ously with radiation produced by two lasers with rela- 
tively broad spectral profiles, which were necessary to 
excite coherently magnetic sublevels that are split by an 
external electric field S (see Fig. [T]) . The model assumes 
that the atoms move classically and are excited at the 
internal transitions. Thus, the internal atomic dynamics 
can be described by a semiclassical atomic density ma- 
trix p, which also depends on the classical coordinates of 
the atomic center of mass. 

The ground state of the Cs atom consists of two hy- 
perfine levels with total angular momentum Fg — i and 
Fg = 4, each containing 2Fg + 1 magnetic sublevels. The 
first laser excites the atoms from the ground state to the 
6P3/2 state, which contains hyperfine levels F^ = 2,3,4, 
and 5. The second laser excites the atom from the 6P3/2 
state to the state, which contained hyperfine levels 

Ff = 1,2,3,4,5, and 6. 

The external electric field £ partially decouples the 
electronic angular momentum from the nuclear spin, 
which implies that the magnetic sublevel energies no 
longer depend quadratically on the electric field (see 
Fig. [T|). In order to obtain the real dependence on the 
electric field, it is necessary to diagonalize the full Hamil- 
ton matrix. It is also necessary to take into consider- 
ation that the decoupling of angular momentum from 
nuclear spin alters the dipole transition probabilities be- 
tween magnetic sublevels. 

The entire model is based on the Optical Bloch Equa- 
tions (OBEs) for the density matrix p (see, for example. 



.^dp 



H,p 



ihRp. 



(1) 



The relaxation operator R includes spontaneous emission 
and transit relaxation. We assume that the density of 
atoms is sufficiently low that different velocity groups of 
thermally moving atoms do not interact. The elements of 
the relaxation matrix are given in Q. The Hamiltonian 
H = Hhfs + V includes the hyperfine Hamiltonian and 
the dipole interaction operator V = — d-E (i), where d is 
the electric dipole operator and E(t) is the electric field 
of the exciting radiation. 

The equations can be simplified by assuming that each 
laser excites only the transition to which it is tuned. We 
also apply the rotating wave approximation for multi- 
level systems [l^ to the OBEs. The resulting stochastic 
differential equations can be further simplified by using 
the decorrelation approach [Tt'I . The stochasticity derives 



from the random fluctuations of the laser radiation with 
finite spectral width. This approach assumes that both 
lasers are uncorrelated and that the integration time for 
each measurement is large compared to the character- 
istic phase-fluctuation time of the exciting light source. 
The decorrelation approximation amounts to solving the 
equations of the density matrix elements that correspond 
to optical coherences and taking a formal statistical av- 
erage over the fluctuating phases This procedure 
results in a system of equations that, when solved, yields 
the observed signals. 

From the density matrix of the final state, one can 
obtain the fluorescence intensities of agiven polarization 
along the unit vector e from [l^ 113, [IH : 
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where /q is a constant and d 
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(5i|d-e|/j) is the 



matrix element between the ground and flnal states of 
the dipole operator along a speciflc polarization direction 
e, i.e., the x— or y— direction. 



III. CALCULATION OF SCALAR AND 
TENSOR POLARIZABILITIES 

A. Motivation 

The description of the signals obtained from the exper- 
iment described above depends on two atomic properties 
simultaneously: the hyperfine constant A and the tensor 
polarizability 02 . If one of these constants can be known 
by independent means, the experiment provides a way 
to determine the other. In this section, we describe an 
all-order relativistic many-body calculation of the tensor 
polarizability a2- A reliable theoretical estimate of this 
constant, together with the experimental results of the 
previous section, can be used to estimate the hyperfine 
constant A, which is difficult to calculate theoretically 
and has not been measured to high precision for the 7,9, 
and IOD5/2 states of cesium. 

B. Method 

The scalar ao and tensor 012 polarizabilities of an atomic 
state V are calculated using formulas 
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where D is the dipole operator and the formula for ao 
includes only the valence part of the polarizability. The 
contribution to ao from the ionic core is negligible for the 
present calculation (16 Oq). The sum over n includes the 
"^3/2 J s-iid ni^5/2 states for the calculation of the 
D3/2 polarizabilities in cesium and the ^^3/2, nFc,/2, and 
n_F'7/2 states for the calculation of the polarizabil- 
ities. The sum over the intermediate states n converges 
rather quickly and only the first few terms need to be 
calculated accurately. Therefore, we separate the calcu- 
lation of the polarizabilities into the calculation of the 
main term q:™^™ and the evaluation of the remainder 
Q,taii^ We include the contributions from the following 
states into the main term: 6P, 7P, 8F, 9P, lOF, IIP, 
12P, 4P, 5F, 6F, IF, and 8F to calculate the polar- 
izabilities of all D states considered in this paper. We 
also include the contributions from the 9P states into 
the calculation of the a™^'"(101?). All electric-dipole re- 
duced matrix elements in Eqs. ^ S]) that are needed for 
the calculation of the main term are calculated using the 
relativistic all-order method, which is briefly described 
below. We use experimental energies from [l^ in the 
main term calculations. We note that the polarizabilities 
of the and lOD states are very sensitive to the val- 
ues of the 9Z) — lOP and 1013 — IIP energy differences, 
respectively, since they are small (50 — 100 cm~^). We 
assume that the energies in Ref. [22] are accurate to all 
quoted digits. The remainders a''^'' are small for all sums 
and are calculated in the Dirac-Hartree-Fock (DHF) ap- 
proximation. 

The all-order method used here sums infinite sets of 
many-body perturbation theory terms. We refer the 
reader to Refs. [H, [l^, [11] for a detailed description of 
the approach. Briefly, the wave function of the valence 
electron v is represented as an expansion 
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where is the lowest-order atomic state function, which 
is taken to be the frozen-core Dirac-Hartree-Fock wave 
function of a state v. This lowest-order atomic state 
function can be written as — aj|0c), where jOc) 

represent DHF wave function of a closed core. The in- 
dices m and n designate excited states and indices a and 
h designate core states. The equations for the excitation 
coefficients are solved iteratively until the correlation en- 
ergy converges to an acceptable accuracy. The excita- 
tion coefficients pma, Pmv, Pmnab, and pmnva are used to 
calculate the matrix elements, which can be expressed 
in the framework of the all-order method as linear or 
quadratic functions of the excitation coefficients. The 
electric-dipole matrix elements as well as the hyperfine 
constants are calculated using the same approach. The 
expansion given by Eq. ^ is restricted to single and dou- 



ble (SD) excitations leading to the omission of certain 
fourth- and higher-order terms. 

We use B-splines [26| to generate a complete set of 
DHF basis orbitals for the all-order calculation. Here, 
we use = 70 splines for each angular momentum. 
The basis orbitals are constrained to a cavity of radius 
R = 220 a.u. The size of the cavity is taken to be large 
enough to fit all of the states needed for the calculation 
of the main terms for all of the polarizabilities calculated 
in this work. The calculation of the polarizabilities of the 
9D and lOD states requires such a large cavity since we 
need to be able to properly describe states up to 12P and 
9P. This work required extensive study of the numerical 
accuracy and stability of the calculations. We verified 
that our basis set gives correct lowest-order (DHF) val- 
ues for the energies of all relevant states and transition 
matrix elements between these states. We have also ver- 
ified that our basis set correctly reproduces DHF values 
of the hyperfine constants for all the nDj states consid- 
ered here. We find that it is necessary to use 70 splines 
to produce an accurate basis set. We also conducted an 
all-order calculation with a smaller cavity {R = 90 a.u.) 
that is appropriate for the calculation of the properties 
of the low-lying states and found that the properties of 
the low-lying states are accurately described by our large 
R = 220 a.u., = 70 basis set. Therefore, we con- 
clude that numerically accurate results can be obtained 
even for such highly-excited states as 12P with the use 
of large basis sets. 



C. Results 

The contributions to the scalar and tensor polarizabil- 
ities for the 71)3/2 state in cesium are listed in Table ID 
We note that the calculation of the scalar and tensor 
polarizability differs only in the angular factor, and all 
matrix elements and energies are the same. The corre- 
sponding energy differences and the absolute values of the 
lowest-order and final all-order electric-dipole reduced 
matrix elements are also listed. The energy differences 
are given in cm~^. Electric-dipole matrix elements are 
given in atomic units (eao), and polarizabilities are given 
in lO'^ Oq, where ao is the Bohr radius. The difference 
between the lowest-order values and the all-order values 
allows to evaluate the size of the correlation correction. 
The accuracy of our calculation is generally higher when 
the relative size of the correlation correction is smaller. 

The contributions from all terms in a™'''" are listed 
separately to identify the most important terms. The 
remainder a*'''' is separated to a''*''(nPi/2), a***'' ("^^3/2)1 
and a''^''(7T.P5/2) for the study of the convergence of these 
three sums. 

We find that three contributions, from the 8P1/2, 
8P3/2J and 5P5/2 states, are dominant. Another term 
(4P5/2) gives a small but significant contribution to the 
tensor polarizability. Therefore, we conduct a more ac- 
curate calculation of the relevant matrix elements and 
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TABLE I: The contributions to scalar and tensor polarizabilities for the 7-D3/2 state in cesium. The corresponding energy 
differences and the absolute values of the lowest-order (DHF) and final all-order electric-dipole reduced matrix elements are 
also listed. The energy differences are given in cm~^. The electric-dipole matrix elements are given in atomic units (eao), and 
the polarizabilities are given in 10"^ flg, where ao is Bohr radius. 
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IOP3/2 


1.642 


1.142 


2706.1 


0.018 


0.014 




IIP3/2 


0.872 


0.627 


3373.2 


0.004 


0.003 




I2P3/2 


0.572 


0.417 


3816.9 


0.002 


0.001 


a'""(nP3/2) 








0.008 


0.006 




4P5/2 


9.165 


13.027 


-1575.4 


-3.9(1) 


0.79(3) 




5P5/2 


46.603 


43.406 


923.7 


74.6(1.1) 


-14.9(2) 




6P5/2 


9.074 


1.289 


2281.9 


0.027 


-0.005 




7P5/2 


5.484 


1.999 


3100.4 


0.047 


-0.009 




8P5/2 


3.767 


1.695 


3631.1 


0.029 


-0.006 


a'""(nF5/2) 








0.434 


-0.087 


Total 










-66.8(1.6) 


71.2(1.2) 



evaluate their uncertainties. The study of the breakdown 
of the correlation correction demonstrates that the main 
contributions to these transitions come from the terms 
containing only single valence excitation coefficients pmv 
(see Eq. (O). In such cases, it is possible to use a semi- 
empirical scaling procedure such as is described, for ex- 
ample, in Ref. [24| to estimate dominant classes of the 
omitted higher-order corrections. The single excitation 
coefficients pmv are multiplied by the ratio of the ex- 
perimental and theoretical correlation energy, and the 
calculation of the matrix elements is repeated using the 
modified excitation coefficients. The difference between 
the ah initio and scaled SD all-order values for the par- 
ticular matrix element is taken to be its uncertainty. The 
relative uncertainty of the corresponding contribution to 
polarizability is twice the relative uncertainty of the ma- 
trix element. As we noted above, we assume that the 
experimental energies are accurate to all digits quoted 
in Ref. ^2^. The uncertainties of the total polarizabil- 
ity values are obtained by adding the uncertainties of the 
individual terms in quadrature. The uncertainty in all re- 
maining contributions is estimated to be insignificant in 
comparison with the uncertainty of the dominant terms. 
We observe significant cancellations between the domi- 
nant terms for both scalar and tensor polarizabilities of 
the 7I?3/2 state. However, the cancellation is more severe 



for the scalar polarizability, where the contributions from 
8-P1/2 a-nd 5F5/2 states are comparable in size but have 
opposite signs. Therefore, we expect higher accuracy of 
our tensor polarizability calculation in comparison with 
the scalar one. 

The contributions to scalar and tensor polarizabilities 
for the state in cesium are listed in Table HIl The 

table is structured in exactly the same way as the one for 
the 7133/2 state. We find that the contribution from the 
8P3/2 state is clearly dominant and the cancellation is 
much less severe. For the tensor polarizability, the next 
largest term, 5i^7 /2 , is six times as small as the dominant 
term. The accuracy of the matrix elements in the dom- 
inant terms is similar for the 7D^/2 ^^nd 7-D5/2 states. 
Therefore, our calculation of the lD^i2 polarizabilities 
is expected to be more accurate than that of the 7D^/2 
polarizabilities. 

The contributions to scalar and tensor polarizabilities 
for the 9Z)j and lODj states in cesium are listed in Ta- 
bles mil and IIVI respectively. The breakdown of the 
polarizability contributions is similar to that of the ID 
polarizability calculations. We list only the dominant 
contributions separately and group all of the other con- 
tributions together in the rows labeled "Other" . The un- 
certainty is evaluated using the method described above. 
The relative importance of the correlation corrections de- 
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TABLE 11: The contributions to the scalar and tensor polarizabilities for the 7-D5/2 state in cesium. The corresponding energy 
differences and the absolute values of the lowest-order (DHF) and the final all-order electric-dipole reduced matrix elements 
are also listed. The energy differences are given in cm~^. The electric-dipole matrix elements are given in atomic units (eao), 
and the polarizabilities are given in 10"^ ap. 



Contribution 


nlj 


7DHF 


7SD 




ao(7-D5/2) 


«2(7D5/2) 




J./ z, 


2.375 


2.909 


-14336.5 


-0.014 


0.014 


7Pl/2 

1-/ ^ 


6.303 


9.679 


-4122.2 


-0.554 


0.554 




8P1/2 


45.594 


43.210 


-277.1 


-164.3(1.7) 


164.3(1.7) 




9^1/2 


16.835 


10.774 


1613.1 


1.755 


-1.755 




lOPl/2 


4.939 


3.555 


2685.1 


0.115 


-0.115 




1 1 P, 

J^J^J 1/2 


2.623 


1 Q47 










1 P 


1 790 
1. ( ZU 




/ yo.y 


U.Ull 


-U.Uil 


„ tail r> \ 








0/17 

U.U4 ( 


0/17 




4^5/2 


2.444 


3.471 


-1596.4 


-0.184 


-0.210 


5F5/2 


12.464 


11.660 


902.7 


3.67(5) 


4.20(5) 




6F5/2 


2.441 


0.457 


2260.9 


0.002 


0.003 




' -^B/2 


1.472 


0.590 


3079.4 


0.003 


0.003 




8F5/2 


1.011 


0.488 


3610.1 


0.002 


0.002 


«'""(nF5/2) 








0.021 


0.024 




4F7/2 


10.925 


15.292 


-1596.5 


-3.6(1) 


1.28(4) 




5Fy2 


55.737 


52.145 


902.6 


73.5(9) 


-26.2(3) 




6F7/2 


10.926 


2.049 


2260.8 


0.045 


-0.016 




7F7/2 


6.588 


2.643 


3079.3 


0.055 


-0.020 




8F7/2 


4.522 


2.186 


3610.1 


0.032 


-0.012 










0.416 


-0.148 


Total 










-89.0(1.9) 


141.8(1.7) 



TABLE IIL The contributions to scalar and tensor polariz- 
abilities for the 9-D3/2 and 9-D5/2 states in cesium in 10^ a%. 



Contribution 00(91)3/2) 02(91)3/2) 



lOPi/2 -1760(9) 1760(9) 

IOP3/2 -483(2) -386(2) 

6P5/2 -129(2) 25.8(4) 

7F5/2 938(8) -188(2) 

Other 31 -22 

Total -1403(12) 1190(10) 



Contribution 00(91)5/2) 02(91)5/2) 



IOP3/2 -2653(12) 2653(12) 

7P5/2 46.3(3) 53.0(4) 

6P7/2 -117(2) 41.9(6) 

7F7/2 927(6) -331(2) 

Other 20 -30 

Total -1777(14) 2386(13) 



TABLE IV: Contributions to the scalar and tensor polariz- 
abilities for the IOD3/2 and IOD5/2 states in cesium in 10'^ a\. 



Contribution 00(101)3/2) 02(101)3/2) 



llPi/2 -4995(24) 4995(24) 

IIP3/2 -1379(6) -1103(5) 

7P5/2 -425(2) 85.1(4) 

8P5/2 2478(16) -496(3) 

Other 84 -65 

Total -4236(29) 3416(24) 



Contribution 00(101)5/2) 02(101)5/2) 



IIP3/2 -7553(31) 7553(31) 

8P5/2 122(1) 140(1) 

7P7/2 -386(3) 138(1) 

8P7/2 2450(17) -875(6) 

Other 51 -89 

Total -5316(36) 6867(32) 



creases with the principal quantum number n and the 
cancellation of different terms becomes less significant 
resulting in smaller uncertainties of the polarizabilities 
for the 913 and 1013 states in comparison with the un- 
certainties for the ID states. Overall, the uncertainties 
of our polarizability calculation are 0.5% — 2.3%. 



D. Comparison with existing experimental values 
and other theory 

Our results for the scalar polarizabilities of the 7-Dj, 
9Z)j, and lODj states in cesium are compared with the 
experimental values from Refs. [lol . [Til . and theoret- 
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TABLE V: Comparison of the scalar polarizabilities ao for 
the 7D, 9D, and lOD states in cesium with other theory and 
experiment. The polarizabihties are given in 10^ ao. 



State 


This work 


Expt. 


Ref. ^3] 


71)3/2 


-66.8(1.6) 


-60(8) [12] 


-65.2 


91)3/2 


-1403(12) 


-1450(120) [llj 


-1400 


IOD3/2 


-4236(29) 


-4185(4) [10] 


-4220 


7D5/2 


-89.0(1.9) 


-76(8) [12] 


-87.1 


9D5/2 


-1777(14) 


-2050(100) [11] 


-1770 


IOD5/2 


-5316(36) 


-5303(8) [10] 


-5300 



TABLE VI: Comparison of the tensor polarizabihties 02 for 
the 7D, 9D, and lOD states in cesium with other theory and 
experiment. The polarizabilities arc given in 10^ ajj. 



State 


This work 


Expt. 


Ref. [13] 


7D3/2 


71.2(1.2) 


74.5(2.0) [4] 


70.4 






66(3) [12J 




9D3/2 


1190(10) 


1183(35) ]4] 


1190 






1258(60) [llj 




IOD3/2 


3416(24) 


3401(4) [10] 


3410 


7D5/2 


141.8(1.7) 


129(4) [12] 


140 


91)5/2 


2386(13) 


2650(140) [llj 


2380 


IOD5/2 


6867(32) 


6815(20) [10] 


6850 



TABLE VII: The breakdown of the correlation correction to 
the hyperfine constants A for the nD-j,/2 and nD^i^ states in 
cesium calculated using the SD all-order method. The ex- 
pressions for all terms are given in [S^]. The values of the 
contributions for the dominant terms and total correlation 
correction are given in % relative to the lowest-order value 
for each state. The total contains contributions from all terms 
{a — t). The normalization factor is also listed. 



Contribution 


5D3/2 


6-D3/2 


7£>3/2 


8D3/2 


9D3/2 


IOD3/2 


Tprm 


11% 


26% 


28% 


28% 


28% 


28% 


Term c 


127% 


57% 


36% 


28% 


23% 


21% 


Term d 


41% 


9% 


4% 


2% 


2% 


1% 


Term h 


13% 


9% 


5% 


4% 


3% 


3% 


Term p 


19% 


13% 


11% 


10% 


9% 


9% 


Total 


214% 


118% 


87% 


75% 


69% 


65% 


Norm 


1.10 


1.14 


1.12 


1.10 


1.10 


1.09 


Contribution 


5D5/2 


6D5/2 


7£>5/2 


8D5/2 


9D5/2 


IOD5/2 


Term a 


-352% 


-264% 


-228% 


-213% 


-205% 


-200% 


Term c 


120% 


54% 


35% 


27% 


23% 


21% 


Term d 


37% 


8% 


3% 


2% 


2% 


1% 


Term h 


-154% 


-28% 


-5% 


3% 


6% 


8% 


Term n 


18% 


16% 


14% 


13% 


12% 


12% 


Term p 


13% 


10% 


9% 


8% 


8% 


8% 


Term r 


-24% 


-18% 


-15% 


-14% 


-14% 


-13% 


Total 


-339% 


-217% 


-184% 


-171% 


-164% 


-160% 


Norm 


1.09 


1.12 


1.10 


1.09 


1.09 


1.08 



ical values from Ref. [H in Table El The polarizabili- 
ties are given in 10^ Oq. The conversion factor from the 
MHz/(kV/cm)^ units to 10'^ atomic units used in the 
present work is 10~''/i/(47reoaQ) = 4.01878, where h is 
the Planck constant. The present values agree with the 
experimental results for 7D3/2, 91)3/2, and IOD5/2 states 
within the corresponding uncertainties. There is some 
discrepancy with the accurate experimental value for the 
101)3/2 state, but the discrepancy is only 1.5 of our es- 
timated uncertainty. However, our values for the 71)5/2 
and 9D5/2 states disagree significantly with the experi- 
mental values for these states. The calculations for the 
7D5/2, 9D5/2, and IOD5/2 state polarizabilities are very 
similar. Thus, the experimental values for the scalar po- 
larizabilities are not consistent with each other accord- 
ing to our theoretical model. Our calculations confirm 
the value for the IOD5/2 state to high precision, and one 
would have expected similar agreement in the case of the 
7D5/2 and 9D5/2 state. 

The results for the tensor polarizabilities for the 7Dj, 
9Dj, and lODj states in cesium are c omp ared with the 
experimental values from Refs. 0, EEllil and theo- 
retical values from Ref. [13] in Table IVTl The polarizabil- 
ities are also given in 10"^ Aq. The present results for the 
nD^j2 states support the measurements of Refs. [1, 
and disagree with the less precise previous measurements 
pH [1^ . The comparison of the nD^ /2 values with ex- 
periment mirrors the result of the comparison for the 
scalar polarizabilities: the 7D5/2 and 91)5/2 values dif- 
fer significantly from the experiment while the IOD5/2 



value agrees with the precise experiment within the cor- 
responding uncertainties. Our values agree with the cal- 
culation of Ref. [l3j for all states for both scalar and 



tensor polarizabilities. 



IV. CALCULATION OF HYPERFINE 
CONSTANTS 

In this section we evaluate the current knowledge 
about the hyperfine constants of the D3/2 and D5/2 states 
of cesium. We describe a calculation of the hyperfine 
constants for the 5D3/2 — IOD3/2 and 5D5/2 — IOD5/2 
states of ^^"^Cs. Then we compare the results of the cal- 
culation to previously measured values. The calculation 
of the hyperfine constants also makes use of the rela- 
tivistic all-order method and is done in the same way 
as the calculation of the electric-dipole matrix elements 
and with the same set of the excitation coefficients Pmai 
Pmv, Pmnab, and Pmnva (scc Eq. ^ ). The breakdown 
of the correlation correction to the hyperfine constants 



A for nD 



and nDc,/2 states in cesium calculated us- 



3/2 

ing the SD all-order method is given in Table IVIIl The 
expressions for the Terms a, c, d, h, n, and p are given in 
[23| . These terms are linear or quadratic functions of the 
excitation coefficients. The values of the contributions 
of the dominant terms and total correlation correction 
are given in % relative to the lowest-order value for each 
state. The normalization factor is also listed. We find 
that the correlation correction is very large, especially 
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TABLE VIIL The hyperfine constants A (MHz) for the nDg/a 
and states in cesium. The lowest-order, "dressed" 

third-order values, and all-order values are compared with 
previous experiments. The experimental data are taken from 



State 


DHF 


_L 1111 H Wl U-Ci. 


All nvri PT" 


F.xnt [51 

J_J A JJ L . \0\ 


5-D3/2 


18.2 


47.0 


52.3 


48.78(7) 




9.27 


21.5 


17.8 


1 6 30^ 1 5") 


7D3/2 


4.70 


10.1 


7.88 


7.4(2) 


8^*3/2 


2.65 


5.46 


4.20 


3.94(8) 


91^3/2 


1.63 


3.28 


2.51 


2.35(4) 


IOD3/2 


1.07 


2.12 


1.62 


1.51(2) 


51)5/2 


7.47 


-32.3 


-16.4 


-21.24(8) 


6D5/2 


3.73 


-8.15 


-3.89 


-3.6(10) 


71)5/2 


1.88 


-2.67 


-1.42 


-1.7(2) 


8D5/2 


1.06 


-1.15 


-0.684 


-0.85(20) 


9D5/2 


0.651 


-0.592 


-0.384 


-0.45(10) 


IOD5/2 


0.428 


-0.343 


-0.238 


-0.35(10) 



for the -D5/2 states where it is several times as large as 
the lowest-order value and has an opposite sign. Ow- 
ing to such an enormous correlation correction, we do 
not expect our results to be very accurate for the 
states. The scaling procedure described above or partial 
ab initio inclusion of the triple excitation as described 
in Ref. [1^ can only evaluate corrections to terms c and 
d, that are not dominant for any of the states except 
51)3/2 ■ Therefore, we can not make an accurate estimate 
of the uncertainty of our values that is independent from 
experimental measurements. 

Our results for the hyperfine constants A (MHz) for the 
nD state in Cs are compared with previous experiments 
in Table IVIIII We list the lowest-order and "dressed" 
third-order values together with the SD all-order values. 
The "dressed" third-order calculation has all lowest-order 
matrix elements replaced by "dressed" matrix elements 
calculated in the random-phase approximation (RPA) 
[2?! . We find large discrepancies between the third-order 
and all-order results indicating very large contributions 
from the fourth- and higher-order terms. Taking into 
account the very large size of the correlation correction 
and obviously large contributions from higher orders, we 
find that the agreement of the all-order calculation with 
measured values is remarkably good. 

We have investigated the issue of the consistency of 
the experimental hyperfine data using our calculation. 
Table IVTIl demonstrates that the breakdown of the corre- 
lation for the 6D — 10D states is rather similar, especially 
for states. We note that nD^j2 and states 

have to be considered separately. The distributions of the 
correlation for both 51)3/2 and 5D5/2 states are clearly 
very different from the ones for the other nD states, and 
these states are omitted from the consistency check be- 
low. For the ^1)5/2 states, the relative contribution of 
Term h changes sign; however, the contribution from this 
term is small in comparison with the experimental un- 



certainty. To cross-check the experimental data, we take 
the experimental value for one particular nDj state and 
rescale it for all the other states with the same J us- 
ing the theoretical values of the correlation corrections. 
The correlation correction is calculated as the difference 
between the final (experimental or theoretical) number 
and the lowest-order DHF value. For example, we take 
the experimental value ^eDg*^ ^^'^ determine how much 
we need to scale our theoretical correlation correction for 
the 61)3/2 state to obtain this value. The scaling factor 
is defined as 



^(6i?3/2 



/jExp- 



3/2 



/I DHF 

^6D3/2 



4SD 



4 DHF ■ 



where ^^hf ^^^^ ^sd lowest-order and all-order 

values from Table IVini for the 6D3/2 state. Next, we 
take our theoretical value for another state, for example. 



ID 



3/2' 



and rescale its correlation correction contribution 



DHF 

ID.,. 



A 



DHF 

ID. in 



using the scaling factor £'(61)3/2): 

A{3D^I2) - 5(6^3/2) X (^fg,^^ - A 

(6) 

Then, we calculate ^(81)3/2), ^(91)3/2), and ^(101)3/2) 
using Eq. ^ . Wc list these values in the column labeled 
5(61)3/2) of Table HXl which indicates that these values 
were obtained with the scaling factor S'(6Z)3/2)- We re- 
peat the procedure using other nD^/2 values to define 
the scaling factor. The uncertainty of the rescaled val- 
ues comes only from the experimental uncertainty of the 
initial experimental value A^-^^^ ■ We find that all results 
in each row are consistent within the uncertainties, lead- 
ing to the conclusion that the experimental results are 
internally consistent. We note that such a procedure will 
not be able to detect a systematic shift of all the experi- 
mental results. Since we cannot accurately evaluate the 
uncertainty of the scaling procedure itself, it is unclear 
if it can yield data that are more accurate than the cor- 
responding experimental data, even though some of the 
rescaled data has smaller uncertainties than the actual 
experimental data. The accuracy of the rescaling is ex- 
pected to be higher when An between the original and 
scaled state is the smallest. 



V. ANALYSIS OF EXPERIMENTAL DATA AND 
ESTIMATE OF THE HYPERFINE CONSTANTS 

The theoretical calculations of the hyperfine constants 
described in the previous section as well as the exper- 
imental measurements of ^Sj contained large uncertain- 
ties. The scaling procedure seems to indicate that the ex- 
perimental values of the review Q, although taken from 
different sources, are consistent with each other. Thus, 
there is an indication that the scaling procedure could 
yield slightly more accurate predictions of hyperfine con- 
stants of states in adjacent levels if the hyperfine constant 
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TABLE IX: The consistency check of the experimental hyperfine constants A (MHz) values for the 710^/2 and states in 

cesium. The actual experimental data from Ref. ^ are listed in the second column. The columns labeled "S(n_Dj)" , n = 6 — 10, 
give data obtained by taking the experimental value for this particular nDj state and rescaling it for all the other states using 
the theoretical values for the correlation correction as explained in the text. The uncertainty of the rescaled values comes only 
from the experimental uncertainty of the initial experimental value nDj. 



State 


Expt. f5] 


S(6I53/2) 


S(7I53/2) 


S(8I53/2) 


S(9I53/2) 


8(10153/2) 


6-03/2 


16.30(15) 




16.5(5) 


16.3(4) 


16.2(4) 


16.1(3) 


7D3/2 


7.4(2) 


7.33(6) 




7.35(16) 


7.3(1) 


7.25(12) 


8153/2 


3.94(8) 


3.93(3) 


4.0(1) 




3.92(7) 


3.89(6) 


9153/2 


2.35(4) 


2.36(2) 


2.38(6) 


2.36(5) 




2.33(3) 


10D3/2 


1.51(2) 


1.53(1) 


1.54(4) 


1.53(3) 


1.52(3) 




State 


Expt. 


S(6i55/2) 


8(7/55/2) 


8(8/55/2) 


S(9I55/2) 


8(10155/2) 


6A5/2 


-3.6(10) 




-4.5(5) 


-4.6(7) 


-4.4(7) 


-5.2(1.1) 


7155/2 


-1.7(2) 


-1.3(4) 




-1.7(4) 


-1.6(3) 


-2.0(5) 


8I55/2 


-0.85(20) 


-0.6(2) 


-0.83(10) 




-0.80(17) 


-0.97(26) 


9155/2 


-0.45(10) 


-0.34(14) 


-0.47(6) 


-0.48(11) 




-0.56(16) 


10155/2 


-0.35(10) 


-0.21(9) 


-0.29(4) 


-0.30(8) 


-0.28(6) 





of one state is known. The experiment described in Sec- 
tion II could provide an independent cross-check of these 
findings. 

With the tensor polarizabilities calculated in section 
III, the simulations described in section II can be used 
to estimate the hyperfine constant. First, we calculate 
a series of simulated curves, varying those experimental 
parameters that we cannot measure precisely, in partic- 
ular the detuning of the lasers. When the overall shape 
of the simulated curve matches the experiment, the po- 
sitions of the features depend on the values of the tensor 
polarizability a2 and the hyperfine constant A. 

We assume that the tensor polarizabilities calculated 
in section III for the 7,9, and IOD5/2 states of cesium 
are the most accurate values available because of the ex- 
cellent agreement between the calculated and previously 
measured values for the IOD3/2 state of cesium. By fix- 
ing the tensor polarizability at the calculated value in our 
simulations, we can thus estimate the hyperfine constant 
A from the level-crossing signals in Figures [H [5l and [6] 

Table |X] summarizes the polarizabilities used in the 
simulations and the hyperfine values A obtained after a 
fit to the experimental data. 

Considering the difficulty in calculating the hfs con- 
stants, the results of the relativistic many-body calcula- 
tion for the hyperfine constant A agree reasonably well 
with the experimental measurements for the 7D5/2 and 
9D5/2 states (within ~ 1.5cr). The large discrepancy in 
the case of the IOD5/2 state seems problematic, since 
the calculations should be internally consistent, if not 
completely reliable in absolute terms. This inconsistency 
could indicate that we slightly underestimated our un- 
certainties. It is also possible that the self-consistency 
check is less reliable in the case of the nD5/2 states, be- 
cause the DHF term and the all order term differ even in 
their sign. 



VI. CONCLUSION 

We obtained new values for the hfs constants A of the 
7,9, and IOD5/2 states. Our values agreed with previously 
measured values, but achieved greater precision. The val- 
ues were obtained by means of measured level-crossing 
signals, a detailed theoretical description of these signals, 
and values for the tensor polarizability calculated with 
an all-order relativistic many-body method. We demon- 
strated the all-order relativistic many-body method's re- 
liability even in highly excited states of ^^'^Cs by com- 
paring scalar and tensor polarizabilities obtained by this 
method with previously experimentally measured values 
for the 7,9,10D3/2 and 7,9,10D5/2 states of ^^^Cs. 

Our calculated polarizability values were in good 
agreement with experiment except for the 7 and 9D5/2 
states. However, the experimental values reported for 
these states are called into question by the fact that 
values reported in the same works for the 7D3/2 [l3 | 
and 9D3/2 [ill] states also disagree with our calculations, 
whereas more recent measurements of the 7 and 9D3/2 
states Q support our calculation, as well as previous cal- 
culations [1J| . The method was further applied to calcu- 
late values for the hyperfine constants A in the 5-IOD3/2 
and 5-IOD5/2 states. Although the calculation cannot be 
considered reliable in absolute terms, nevertheless they 
agreed reasonably well in the case of the 7D5 /2 and 9D5 /2 
states. For the IOD5/2 state, the agreement was not as 
good. 
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TABLE X: Comparison of the experimentally obtained hyperfine constants with previous experiments and the theory presented 
in this work 



Cesium 
atomic 
state 


Calulated 
tensor 
polarizability 
(lO^ag) 




hyperfine constant 
(MHz) 




This work 


Previous 
experiment 


Theory 


7D5/2 


141.8(1.7) 


-1.56(9) 


-1.7(2) [5] 


-1.42 


9D5/2 


2386(13) 


-0.43(4) 


-0.45(10) [5] 


-0.384 


IOD5/2 


6867(32) 


-0.34(3) 


-0.35(10) [5] 


-0.238 
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